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Abstract 

The potential of a linear e + e~ collider operated at a centre-of-mass 
energy of 350 GeV is studied for the measurement of the Higgs boson 
mass. An integrated luminosity of 500 fb" 1 is assumed. For Higgs 
boson masses of 120, 150 and 180 GeV the uncertainty on the Higgs 
boson mass measurement is estimated to be 40, 65 and 70 MeV, re- 
spectively. The effects of beam related systematics, namely a bias 
in the beam energy measurement, the beam energy spread and the 
luminosity spectrum due to beamstrahlung, on the precision of the 
Higgs boson mass measurement are investigated. In order to keep 
the systematic uncertainty on the Higgs boson mass well below the 
level of the statistical error, the beam energy measurement must be 
controlled with a relative precision better than 10~ 4 . 



1 Introduction 

In the Standard Model [1] particles acquire mass due to spontaneous sym- 
metry breaking by introducing a doublet of complex scalar fields. This so 
called Higgs mechanism [2] leads to one scalar particle, the Higgs boson. 
The mass of the Higgs boson is a free parameter of the Standard Model 
and of fundamental nature. If the Higgs boson exists, the Large Hadron 
Collider at CERN will be able to discover it [3]. Precision measurements 
of the Higgs boson parameters and the exploration of the complete Higgs 
boson profile will be one of the central tasks at a future linear e + e~ collider. 

In this article we study the potential of a future e + e~ collider for the 
measurement of the mass of a relatively light Higgs boson, in the mass 
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range from 120 to 180 GeV, and investigate possible systematic effects in- 
fluencing the precision of this measurement. The analysis presented ex- 
tends previous studies on the measurement of the mass of a light Higgs 
boson [4] and complements recent studies on the determination of reso- 
nance parameters of the Higgs boson with the mass in the range from 200 
to 320 GeV at a future linear e + e~ collider [5]. 



2 Experimental Conditions and Detector Simu- 
lations 

The study is performed for a linear collider operated at a centre-of-mass 
energy, s/s, of 350 GeV and an event sample corresponding to an inte- 
grated luminosity of 500 fb _1 . This integrated luminosity is expected in 
about one year of running with design luminosity 

The detector used in the simulation follows the proposal for the TESLA 
collider presented in the Technical Design Report [6]. The interaction re- 
gion is surrounded by a central tracker consisting of a silicon micro-vertex 
detector as the innermost part and a time projection chamber. In radial 
direction follow an electromagnetic calorimeter, a hadron calorimeter, the 
coils of a superconducting magnet and an instrumented iron flux return 
yoke. The solenoidal magnetic field is 4T. The central tracker momentum 
resolution is 

^ = 7-KT 5 - Pt , (1) 

Pt 

where p t is the transverse momentum in GeV/ c. 

The energy resolutions of the electromagnetic and hadron calorimeters 
are: 

^=i§®0.6%, ^ = ^®4%, (2) 

where E e and Eh are the energies of electrons and hadrons in GeV. The 
polar angular coverage of the central tracker maintaining the resolution is 
| cos^l < 0.85, above this range the tracking resolution deteriorates. The 
electromagnetic and hadron calorimeters cover | cos#| < 0.996 maintain- 
ing the resolution over the whole angular range. The simulation of the 
detector is done using the SIMDET [7] package. 

The event reconstruction is done in terms of particle flow objects. First, 
tracks are measured with the tracking system and associated to calorime- 
ter clusters to define charged particle flow objects of electrons, muons and 
charged hadrons. Since the momentum measurement by the tracking sys- 
tem is much more accurate than the angular and energy measurements 
with calorimeters, the tracking information is used for the determination 
of the four-momentum of charged particles. Calorimetric clusters with no 
associated track are regarded as neutral particle flow objects originating 
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from photons and neutral hadrons. Measurements of the four-momentum 
of neutral objects are solely based on the calorimetric information. 

3 Physics Processes 

At a centre-of-mass energy of ^fs = 350 GeV, the dominant process for light 
Higgs boson production in the Standard Model is e + e~ — ► ZH. Events of 
this process, hereafter referred to as signal, are generated using PYTHIA 
[8] for Higgs boson masses, m H , of 120 GeV, 150 GeV and 180 GeV. For the 
Higgs boson, all decay modes are simulated as expected in the Standard 
Model. The decay modes into hadrons, WW and ZZ are investigated in 
detail. Z decays are considered into electrons, muons and hadrons. The 
Standard Model cross sections are given in Table [l] for the investigated 
signal channels. 



decay mode 


topology 


cross section, fb 






120 GeV 


150 GeV 


180 GeV 


ZH -> £+£-qq, ZH -> £+£"gg 


2£+2-jets 


8.8 


2.0 


0.06 


ZH — > qqq'q', ZH — > qqgg 


4-jets 


91.9 


20.4 


0.62 


ZH -> f+rWW, W -> qq' 


2£+4-jets 


0.6 


2.6 


2.6 


ZH -> qqWW, W -> qq' 


6-jets 


6.0 


26.5 


26.6 


ZH -> £ + £~ZZ, Z -> qq 


2£+4-jets 


0.08 


0.33 


0.17 


ZH — > qqZZ, Z — > qq 


6-jets 


0.82 


3.46 


1.73 


ZH — > all final states 




160.3 


123.7 


89.0 



Table 1: The cross sections, in fb, times the branching fractions of the inves- 
tigated signal final states for Higgs boson masses of 120 GeV, 150 GeV and 
180 GeV as predicted in the Standard Model for s/s = 350 GeV. Also given 
is the total cross section for e + e~ — > ZH. The cross sections are calculated 
with PYTHIA taking into account initial state radiation. 



For background estimations events are generated with PYTHIA for the 
processes e + e~ — > qq(7), e + e" — > W + W~, e + e~ — > Z(7*)Z(7*) and e + e~ — > 
7*7*e + e~ — > ffe + e~. Six fermion final states resulting from the triple gauge 
boson production are generated with the WHIZARD package [9]. The 
cross sections of the main background reactions are given in Table |2j The 
numbers of events generated for each background channel as well as the 
number of generated signal events correspond to an integrated luminosity 
of 500 fb 1 . 

Initial state radiation is simulated by PYTHIA. Beamstrahlung is taken 
into account using the CIRCE program [10]. 
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background process 


cross section, fb 


events 


e + e" — > 7*7*e + e" — > ffe + e" 


4.0 x 10 6 


2.0 x 10 9 


e + e~ — > qq(7) 


2.7 x 10 4 


1.4 x 10 7 


e+e" -> W+W" 


1.3 x 10 4 


6.5 x 10 6 


e+e" -> Z(7*)Z( 7 *) 


1.0 x 10 3 


5.0 x 10 5 


e+e" -> W+W-Z 


13.2 


6.5 x 10 3 


e+e" -> ZZZ 


0.48 


9.6 x 10 2 



Table 2: The cross sections, in fb, and the numbers of events expected for 
the important background processes at a centre-of-mass energy 350 GeV. 
An integrated luminosity of 500 fb -1 is assumed. 



4 Analysis Procedure 

The measurement of the Higgs boson mass is based on the reconstruction 
of the exclusive final states ZH -> £ + £~qq, ZH -> qqq'q', ZH -> £ + £~WW 
and ZH — > qqWW. In the latter two cases, the contributions from the ZH — > 
£ + £~ZZ and ZH — > qqZZ final states are also taken into account. 

The analysis in all channels proceeds as follows. First, a selection of 
events of the specific topology is applied to the samples of signal and back- 
ground events exploiting event shape variables and lepton identification. 
For Higgs boson masses below 150 GeV, the decay H — > bb is dominant, 
leading to 2£+2-jet and 4-jet topologies. For larger Higgs boson masses, 
the H — > WW decay becomes dominant, leading to 2^+4-jet and 6-jet final 
states. For their study, jet identification is crucial. A kinematic fit imposing 
energy and momentum conservation improves considerably the di-jet or 
4-jet mass resolutions, and hence the accuracy of the mass measurement. 
In order to construct the covariance matrix used in the kinematic fit, the 
resolutions of the lepton and jet energies and angular measurements are 
needed. For the lepton momentum, the tracker resolution from Equation 
(1) is used. The resolution in the polar and azimuthal angles, 9 and 0, of 
the lepton momentum vector are obtained from Monte Carlo as: 

oq = 1 mrad, = (3) 

The resolutions of the jet energies and angular measurements are obtained 
from a Monte Carlo study using the sub-detector resolutions from Equa- 
tions (1) and (2). They are parameterized as: 

IT? 30% it; A ae (AS 

a E E = —=, 00 = 15 mrad, as = -— -, (4) 

y/E smO 



where E is the energy, 9 the polar and the azimuthal angle of the jet. 

The Higgs boson mass is obtained by fitting the invariant mass spec- 
trum of the jets assigned to the Higgs boson decay with a superposition of 
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signal and background distributions. The shape of the signal distribution 
is derived from a high statistics Monte Carlo sample of signal events and 
kept fix in the fit. Free parameters are the peak value and the normalisa- 
tion factor of the signal distribution. 

4.1 The ZH E+l'qq and ZH -» qqq'q' Final States 

These final states are characterized by two isolated leptons and two jets 
or by four jets and have the full energy deposited in the detector. Hence, 
events where the total energy visible in the detector is less that 80% of the 
centre-of-mass energy are rejected. Global event characteristics are used 
for the signal selection. 

For the channel ZH — > £ + £~qq, the number of reconstructed particles 
must be greater than 20, the event thrust, T, must be less than 0.85 and the 
absolute value of the cosine of the polar angle of the thrust vector, cos 9t, 
must be less than 0.9. Electrons are identified as energy deposits in the 
electromagnetic calorimeter whose shape is compatible with the expecta- 
tion for an electromagnetic shower and with a matched track in the central 
tracker. The measured track momentum and shower energy must be in 
agreement within 5% and the shower leakage into the hadron calorimeter 
must be less than 2 GeV. Muons are tracks pointing to energy deposits 
in the calorimeters which are consistent with the expectation for a mini- 
mum ionizing particle. A pair of electrons or muons with opposite charge 
is required. Both electrons and muons must have momenta larger than 
10 GeV/ c and fulfill the polar angle cut | cos 8 e \ < 0.9. Leptons must satisfy 
isolation criteria, meaning that there are no other particles reconstructed 
within a 15° cone with respect to the lepton momentum vector. The in- 
variant mass of a pair of leptons must be compatible with the mass of the 
Z boson within 10 GeV. These criteria reduce the backgrounds listed in 
Table in the selected sample to the level of a few % with the exception of 
the process e + e _ — >• ZZ. A cut on the polar angle of the momentum vector 
of the di-electron or di-muon system, | cos 9 U \ < 0.9, further suppresses the 
ZZ background. The signal selection efficiency is about 45%. All recon- 
structed particles, except the two isolated leptons, are grouped into two 
jets using the Durham [11] jet clustering algorithm. 

Event selection for the ZH — > qqq'q' channel is performed by requiring 
the number of reconstructed particles to be larger than 40, T < 0.85 and 
| cos8 T \ < 0.8. No isolated leptons with an energy greater than 10 GeV 
are allowed. Reconstructed particles are grouped into four jets using the 
Durham jet clustering algorithm. Events are retained if the jet resolution 
parameter, for which the event is resolved from the four- to three-jet topol- 
ogy, yu, fulfill the relation log (y 34 ) > —5. 

The selected events of both final states are subject to a kinematic fit [12] 
imposing energy and momentum conservation. The kinematic fit is per- 
formed by varying the lepton momenta and angles within their resolu- 
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tions given by Equations (1) and (3), respectively. The jet energies and an- 
gles are varied within the corresponding resolutions given by Equations 
(4). 

For events selected as ZH — > £ + i qq, energy and momentum conserva- 
tion results in four constraints (4C fit). Since the experimental resolution 
in the invariant mass of the di-lepton system is much smaller than the nat- 
ural width of the Z boson, no constraint is applied in the kinematic fit to 
force the di-lepton mass to rriz- The di-jet invariant mass spectra after the 
4C fit are shown in Figure |l] for ran = 120 GeV and 150 GeV, respectively. 
Clear signals are seen on top of the remaining smooth background from 
e + e~ — > ZZ. Also shown are the contributions from H — > WW and H — > LL 
decays to the signal. These are negligible for m H = 120 GeV but amount to 
62% and 5%, respectively of the signal for m H = 150 GeV. 

The masses obtained from the fits equal the generated Higgs boson 
masses and have errors of 85 MeV for m H = 120 GeV and 100 MeV for m H 
= 150 GeV. 

For the 4-jet final states, in addition to the four constraints from energy 
and momentum conservation, the invariant mass of the two jets assigned 
to the Z boson decay is constrained to m z . Hence, a 5C fit is performed 
for all possible di-jet pairings. The pairing with the minimal x 2 is chosen. 
In addition, this x 2 must be less than 70. The signal selection efficiency 
is about 25%, however the remaining event sample contains consider- 
able background from e + e~ — > TIL, e + e~ — > W + W~ and e + e~ — > qq(7). The 
signal-to-background ratio is enhanced using the identification of b-quark 
jets. The ZVTOP [13] topological vertex finder adapted for the pixel micro- 
vertex detector [6] is used to search for secondary vertices inside jets and 
determine mass, momentum and decay length of the vertex. In addition, 
the impact parameter joint probability [14] and the two highest impact pa- 
rameter significances are used as input into neural networks trained with 
jets containing no, one and more than one secondary vertices. A jet b-tag 
variable is defined [15] as function of the neural network output x as 



fb{x) + fudscW 

where fb and f u d sc are probability density functions of the neural network 
outputs in samples of b-jets and udsc-jets, respectively. The improvement 
in the signal to background ratio can be seen in Figure |2j where invariant 
mass distributions are shown for the two jets assigned to the Higgs bo- 
son decay without a requirement on the jet b-tag variable and for events 
with at least two jets satisfying B(x) > 0.2. The latter requirement keeps 
the signal statistics almost unchanged but reduces significantly the back- 
ground. For mu = 150 GeV also events from ZH — > qqZZ and ZH — > qqWW 
are selected amounting to 25% of the signal. Using the invariant mass dis- 
tributions obtained with the requirement B(x) > 0.2 the fit of the Higgs 
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boson mass is performed. As an example, the di-jet invariant mass distri- 
bution and the fitted function of the signal is shown in Figure |3j 

The results for the Higgs boson masses are equal to the generated masses. 
The statistical errors are 45 MeV at mu =120 GeV and 170 MeV at mu = 
150 GeV. 

4.2 The ZH £+£~WW and ZH -»> qqWW Final States 

We consider W-boson decays into two quarks, hence the topologies of 
these final states are two isolated leptons accompanied by four jets or six 
jets, respectively. The requirements for electron and muon identification 
are the same as in the previous section. Although event selection is opti- 
mized specifically for the ZH — > £ + £~WW and ZH — > qqWW final states, 
contributions from the ZH — > £ + £~ZZ and ZH — ► qqZZ channels are also 
taken into account. 

Events are selected with an energy deposited in the detector of more 
than 80% of the centre-of-mass energy and a number of the reconstructed 
particles larger than 40. 

Events of the final state ZH — > £ + £~WW must contain a pair of isolated 
electrons or muons with opposite charges. Furthermore, the event thrust 
and the polar angle of the thrust vector are used to suppress the dominant 
background from the WW and ZZ final states. The values of the cuts are 
T < 0.95 and | cos# T | < 0.95. Since the two leptons of ZH -> £ + £~WN 
originate from the Z decay, their invariant mass is required to be equal 
within 10 GeV to mz- A cut on the polar angle of the di-lepton momentum 
vector, | cos 0«| < 0.9, further suppresses the ZZ background. Tracks and 
calorimetric energy deposits not stemming from the leptons are grouped 
into four jets using the Durham algorithm. The jet resolution parameter 
I/34 must satisfy log(y 34 ) > —6.0. 

Then a 4C kinematic fit is performed imposing energy and momentum 
conservation. Only events for which the \ 2 of the 4C fit is less than 50 are 
retained in the selected sample. The signal selection efficiency amounts 
to 50% at m H = 150 GeV and 60% at m H = 180 GeV. The 4-jet invariant 
mass distributions after the kinematic fit are shown in Figure |U for m H = 
150 GeV and 180 GeV. 

The fit of the mass spectra in Figure HI again results in mass values for 
the Higgs boson equal to the generated ones. The uncertainties of the 
masses amount to 90 and 80 MeV for m H = 150 GeV and m H = 180 GeV, 
respectively. 

The small background in this channel comes mainly from the semilep- 
tonic decays of pair produced Z bosons and triple gauge boson produc- 
tion, ZWW, with a leptonic Z decay. Events of the process ZH — » £ + £~ TIL 
constitute 13% and 6% of the signal in the selected sample for mu = 150 
and 180 GeV, respectively. 

The ZH — > qqWW channel is selected by requiring T < 0.9 and | cos 9 T \ 
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< 0.95. There must be no isolated leptons with an energy greater than 
10 GeV. The reconstructed particles are grouped into six jets using Durham 
jet algorithm. The jet resolution parameter, for which an event is resolved 
from the 6- to 5-jet topology, y^, must satisfy \og(y 5e ) > —8. Then a likeli- 
hood discriminant, L H z/ is defined using as input the number of particles 
reconstructed in an event, the polar angle of the thrust vector and the jet 
resolution parameters y 34 and y 56 . Events are accepted when the value of 
this discriminant is larger than 0.9. As an example, Figure [5]shows the dis- 
tribution of Lhz for the signal events for mu = 180 GeV and the background 
processes. The six jets are now grouped in three di-jet pairs following crite- 
ria which depend on the mass of the Higgs boson. For mu < 2m w usually 
only one W is expected to be on the mass shell, while the other is produced 
with a mass close to the difference between m H and m w . The quantity 

X 2 = (my - m z ) 2 /al + (m k i - m w ) 2 /Vw + ( m mn - m klmn + m w ) 2 /o"w* 

is calculated for all possible di-jet combinations, where my is the invariant 
mass of the two jets assigned to the Z boson, m k i the invariant mass of 
two jets assigned to the on-shell W boson, m mn the invariant mass of two 
jets assigned to the off-shell W boson and m klmn the invariant mass of the 
four jets assigned to decay H — > WW*. The quantities a\, and cx^* 
are obtained from Monte Carlo studies as the convolution of the bosonic 
widths and the mass resolutions and are estimated to be 6, 9 and 15 GeV, 
respectively. For mn > 2mw both W bosons are on shell. Hence all di-jet 
combinations are taken and the quantity 

X 2 = (my - m z ) 2 /o"z + ( m ki - mw) 2 /°"w + (m mn - mw) 2 /^w 

is calculated. The jet pairing with the smallest value of x 2 is chosen and 
subject of a kinematic fit imposing energy-momentum conservation and 
constraining the mass of the two jets assigned to the Z boson to m z . Events 
are selected into the final sample if the x 2 °f the 5C fit is less than 30. In 
addition, the fitted mass of the jets originating from the on-shell W decay 
must be equal to m w within 20 GeV in the event sample selected for m H = 
150 GeV. For m H = 180 GeV, the sum and the difference of the fitted masses 
of the two jet pairs assigned to a W decay must be between 125 GeV and 
185 GeV and —20 and 20 GeV, respectively. The signal selection efficiency 
amounts to about 20%. The sample selected for m H = 150 GeV also contains 
5% signal from the ZH — > qqq'q' final state. 

The distribution of the invariant mass of the 4-jet system is shown in 
Figure |S] for m H = 150 GeV and 180 GeV, respectively. From the fit ap- 
proximating the signal by a Gaussian the uncertainties of the masses are 
100 MeV and 150 MeV for m H = 150 GeV and 180 GeV, respectively. The 
background in this channel originates from e + e~ — > W + W~, e + e~ — > TIL 
and e + e~ — > qq(7) final states, and from triple gauge boson production 
processes. Events of the process ZH — > qqZZ constitute 9% and 5% of the 
signal in the selected sample for m H = 150 and 180 GeV, respectively. 
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4.3 Combined Results 



Table |3] summarizes the statistical accuracy on the determination of mn for 
the different final states and their combination. It should be noted that 
considerable overlap exists in the selected samples of the ZH — > £ + £~qq 
and ZH -> £+£~WW channels and of the ZH -> qqq'q' and ZH -> qqWW 
channels. Hence, the combination is performed only for the non-overlapping 
topologies which gives a minimal combined error on the Higgs boson 
mass. This is done using the formula: 

A 2 (m H ) Y A?(m H )' 

where A is the combined error, whereas Aj is the error obtained in the i th 
channel. 



Decay mode 


A(m H ) in MeV 


120 


150 


180 


ZH -> £+£"qq 


85 


100 




ZH — > qqq'q' 


45 


170 




zh -> frww 




90 


80 


ZH -> qqWW 




100 


150 


Combined 


40 


65 


70 



Table 3: Uncertainties on the determination of the Higgs boson mass for 
m H = 120, 150 and 180 GeV. The ZH -> £+£~WW and ZH -> qqWW chan- 
nels are used for the combination at m H = 150 GeV. 



5 Beam Related Systematic Effects 

We have investigated the effect of a bias in the beam energy measurement, 
of the beam energy spread and of an uncertainty in the differential lumi- 
nosity spectrum on the measurement of the Higgs boson mass. 

The impact of a bias in the beam energy measurement is estimated by 
generating signal samples with both positron and electron beam energies 
shifted with respect to the nominal value of \fs/2. These shifts are varied 
from -100 MeV to 100 MeV in 25 MeV steps. Since in the kinematical fit the 
energy is constrained to the nominal value, y/s= 350 GeV, the shift in the 
beam energy is expected to result in a shift in the measured Higgs boson 
mass. 

As an example Figure shows the distributions of fitted values of m H 
in the ZH — > £ + £~qq channel for shifts in the beam energies of +25 MeV, 
zero MeV and —25 MeV. In each of the three considered cases the distribu- 
tion of m H is obtained from 200 statistically independent signal samples. 
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The shift obtained in the fit of mn corresponds roughly to the shift of the 
beam energy with opposite sign. In the range of beam energy shifts from 
-100 to 100 MeV the shift in the Higgs boson mass is found to depend lin- 
early on the shift in the beam energy: 

5m H = —a ■ SE h , (5) 

with a = 0.85 for the ZH -> qqq'q' channel, 0.80 for the ZH -> qqWW chan- 
nel, and 1.04 for the ZH -> £+£~qq and ZH -»■ £+£~WW channels. Hence, 
in order to keep the systematic bias in m H well below its statistical error, 
the beam energy measurement must be controlled with a precision better 
than 10" 4 . 

To estimate the impact of a beam energy spread, a Gaussian distribu- 
tion of the beam energy has been used for the generation of signal events. 
As an example, Figure El shows the reconstructed Higgs boson mass spec- 
trum for a sample of ZH — > qqq'q' events for a 1% energy spread for both 
electron and positron beams and the same distribution for a fix beam en- 
ergy of yfs/2. The inclusion of this beam energy spread slightly broadens 
the mass distribution and degrades the precision obtained for m H from 
45 MeV to 50 MeV in the ZH -> qqq'q' channel and from 85 MeV to 90 MeV 
in the ZH — > £ + £~qq channel. For the TESLA machine the expected energy 
spread amounts to 0.15% for the electron beam and 0.03% for the positron 
beam [6]. For these values of the beam energy spread, no significant degra- 
dation of the precision in the Higgs boson mass measurement is observed. 

The energy spectra of the colliding electrons and positrons at a high 
energy linear collider will be significantly affected by photon radiation 
of the particles in one bunch in the coherent field of the opposite bunch. 
This effect is referred to as beamstrahlung. The program CIRCE provides 
a fast simulation of the beamstrahlung under the assumptions that the 
beamstrahlung in the two beams is equal and uncorrelated between the 
beams. The spectrum is parameterized according to 

f(x) = a 5(l -x) + aix a2 (l - x) a \ 

where x is the ratio between the energies of the colliding electron and 
positron and the initial energy of the undisrupted beam. The parameters 
Oj depend on the operational conditions of the linear collider. The nor- 
malization condition, J f{x)dx = 1, fixes one of these parameters, leaving 
only three of them independent. The default parameters for the TESLA 
machine operated at the centre-of-mass energy of 350 GeV are: 

a = 0.55, ai = 0.59, a 2 = 20.3, a 3 = -0.63. 

It has been shown that from the analysis of the acollinearity spectrum of 
Bhabha scattering events, the parameters can be determined with a pre- 
cision of about 1% [16]. To visualize the effect of the uncertainty in the 
determination of the parameters a if the beam energy spectra are shown 
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in Figure |9]f or nominal values of the parameters a.; and for the parameter 
a shifted by ±10% from its nominal value. Figure llOl presents the corre- 
sponding Higgs boson mass spectra a the sample of ZH — > £ + £~qq events. 
An uncertainty of 10% in the determination of the parameters ao results 
in a systematic uncertainty of about 10 MeV on the Higgs boson mass in 
the ZH — > £ + £~qq and ZH — > qqq'q' channels. The same result is obtained 
for the other parameters. The uncertainty is reduced to about 1 MeV if the 
parameters a« are measured with an accuracy of 1%. The same result is 
obtained for the study of the ZH -> £ + £~WW and ZH -> qqWW channels. 

6 Conclusion 

The potential of the future linear e + e collider for the measurement of the 
Higgs boson mass is evaluated. Assuming an integrated luminosity of 500 
fb _1 , the Higgs boson mass can be measured with a statistical accuracy 
ranging from 40 MeV to 70 MeV for m H between 120 GeV and 180 GeV. In 
order to keep the systematic uncertainty due to a bias of the beam energy 
measurement well below the statistical uncertainty, the beam energy mea- 
surement has to be controlled with a precision better than 10~ 4 . Under 
operational conditions envisaged for the TESLA machine, the beam en- 
ergy spread and uncertainty in the differential luminosity spectrum are 
found to have negligible effect on the Higgs boson mass measurement. 
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Figure 1: The di-jet invariant mass from the ZH — ► £ + £ qq final state after 
a 4C kinematic fit for m H = 120 GeV (top) and 150 GeV (bottom). 
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Figure 2: The distribution of the invariant mass of the two jets assigned to 
the Higgs boson decay in the ZH — > qqq'q' final state without requirement 
on the b-tag and after requiring the values of the b-tag of two jets to be 
larger than 0.2. 
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Figure 3: The invariant mass of the two jets assigned to the Higgs boson 
decay in the ZH — > qqq'q' final state after the 5C kinematic fit for m H = 
120 GeV. For two jets the b-tag must be larger than 0.2. 
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Figure 4: The 4-jet invariant mass from the ZH — > £ + £ WW final state after 
a 4C kinematic fit for m H = 150 GeV (top) and 180 GeV (bottom). 
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Figure 5: The distributions of the signal likelihood used to select 
H — » WW — > 6-jet final states for m H = 180 GeV. Solid, dashed and dotted 
lines represent the background processes, ZH — > qqWW and ZH — > qqZZ 
signals. The vertical line indicates the cut imposed on this quantity. 
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Figure 6: The four jet invariant mass from the ZH — > qqWW final state 
after a 5C kinematic fit for m H = 150 GeV (top) and m H = 180 GeV (bottom). 
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Figure 7: The spectrum of the fitted values of the Higgs boson mass as 
obtained from 200 independent signal samples for the case when both 
electron and positron beam energies are overestimated by 25 MeV (dotted 
histogram), when they are underestimated by 25 MeV (dashed histogram) 
and when no shifts are introduced to the beam energies (solid histogram). 
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Figure 8: Reconstructed Higgs boson mass spectrum in the sample of 
the ZH — > qqq'q' events for the case of monochromatic beams (solid his- 
togram) and for the case of 1% Gaussian energy spread for both electron 
and positron beams (dashed histogram). 
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Figure 9: The beam energy spectrum after beamstrahlung for nominal pa- 
rameters di at y/s = 350 GeV (solid histogram) and for the cases when the 
parameter a is shifted from its nominal value by -10% (dashed histogram) 
and +10% (dotted histogram). E c is the energy of the colliding particle in- 
cluding beamstrahlung and E bea m is the nominal beam energy. 




Figure 10: The reconstructed Higgs boson mass spectrum in the sample 
of ZH — > £ + £~qq events for nominal parameters dj (see text) at 350 GeV 
centre-of-mass energy (solid histogram) and for the cases when parameter 
a is shifted from its nominal value by -10% (dashed histogram) and +10% 
(dotted histogram). 
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